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Abstract 

Ergothioneine is a small, sulfur-containing metabolite (229 Da) synthesized by various species of bacteria and fungi, which 
can accumulate to millimolar levels in tissues or cells (e.g. erythrocytes) of higher eukaryotes. It is commonly marketed as a 
dietary supplement due to its proposed protective and antioxidative functions. In this study we report the genes forming 
the two-step ergothioneine biosynthetic pathway in the fission yeast, Schizosaccharomyces pombe. We identified the first 
gene, egtl^ (SPBCl 604.01), by sequence homology to previously published genes from Neurospora crassa and 
Mycobacterium smegmatis. We showed, using metabolomic analysis, that the Aegtl deletion mutant completely lacked 
ergothioneine and its precursors (trimethyl histidine/hercynine and hercynylcysteine sulfoxide). Since the second step of 
ergothioneine biosynthesis has not been characterized in eukaryotes, we examined four putative homologs (Nfsl/ 
SPBC21D10.11C, SPAC11D3.10, SPCC777.03c, and SPBC660.12c) of the corresponding mycobacterial enzyme EgtE. Among 
deletion mutants of these genes, only one {ASPBC660.12c, designated Aegt2) showed a substantial decrease in 
ergothioneine, accompanied by accumulation of its immediate precursor, hercynylcysteine sulfoxide. Ergothioneine- 
deficient strains exhibited no phenotypic defects during vegetative growth or quiescence. To effectively study the role of 
ergothioneine, we constructed an egtl^ overexpression system by replacing its native promoter with the nmtT^ promoter, 
which is inducible in the absence of thiamine. We employed three versions of the nmtl promoter with increasing strength 
of expression and confirmed corresponding accumulations of ergothioneine. We quantified the intracellular concentration 
of ergothioneine in S. pombe (0.3, 157.4, 41.6, and up to 1606.3 iiM in vegetative, nitrogen-starved, glucose-starved, and 
egf /^-overexpressing cells, respectively) and described its gradual accumulation under long-term quiescence. Finally, we 
demonstrated that the ergothioneine pathway can also synthesize selenoneine, a selenium-containing derivative of 
ergothioneine, when the culture medium is supplemented with selenium. We further found that selenoneine biosynthesis 
involves a novel intermediate compound, hercynylselenocysteine. 
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Introduction 

Ergothioneine (EGT) is a sulfur-containing jVb(jjVo;,jVa-trimethyI- 
L-histidine-derived metabolite that is synthesized by various 
species of bacteria and fungi (recently extensively reviewed by 
Cheah and Halliwell [1]). Higher organisms obtain EGT in food 
and accumulate it in certain tissues up to millimolar levels [2,3] 
through a specific transporter, ETT/OCTNl [4]. In mammals, 
large amounts of EGT are found in erythrocytes, bone marrow, 
liver, kidney, eye lens, and seminal fluid [2,5,6]. Nevertheless, 
EGT is neither a nutrient (it is virtually unmetabolized in humans) 
nor a vitamin (it is non-essential). EGT is commonly marketed as a 
dietary supplement or nutraceutical, due to its anti-oxidant 
properties in vitro, reported in numerous publications [3,7-10]. 
Direct scavenging of free radicals and chelation of transition 



metals are the most widely cited possible functions of EGT 
[10,11]. However, so far, no rigorous research has conclusively 
demonstrated any benefit of EGT in vivo. It is unclear whether 
EGT consumption contributes to human health, and if it does, 
what daily intake is optimal. It is thus of considerable interest for 
biologists and medical scientists to uncover the physiological 
mechanism of EGT at the molecular level. Recently, biosynthetic 
pathways for EGT have been characterized in Mycobacterium 
smegmatis [12] SlXiA Meurospora crassa [13], allowing the use of genetic 
methods. 

The fission yeast, Schizosaccharomyces pombe, is a suitable model 
organism for the study of cell division and quiescence [14-16]. We 
have previously established a method of comprehensive metabo- 
lomic analysis in S. pombe using liquid chromatography-mass 
spectrometry (LC-MS) [17]. Among the several hundred observed 



PLOS ONE I www.plosone.org 



1 



May 2014 | Volume 9 | Issue 5 | e97774 



Ergothioneine in Fission Yeast 



metabolites, we also identified EGT and described its accumula- 
tion under glucose starvation [18] and nitrogen starvation [19]. In 
addition, we reported abnormally high accumulations of EGT in 

the proteasome regulatory subunit mutant mts3-l, which suffers 
fi'om severe oxidative stress caused by mitochondrial dysfunctions 
in GO arrest [20] , and in the Kriippel-like zinc-finger transcription 
factor deletion mutant Aklfl, which exhibits cell wall defects in 
long-term quiescence, accompanied by up-regulation of mito- 
chondrial transcripts [21]. Interestingly, EGT is not synthesized by 
the budding yeast, Sacchammyres cerevisiae [22], suggi^sting that not 
all fungi require EGT for normal ccllukir function. Due to its 
simplicity and advanced genetics, S. pombe might represent an ideal 
unicellular, eukaryotic system to study the biochemical role of 
EGT. 

Here we report identification of genes forming the two-step 
biosynthetic pathway of EGT from histidine in S. pombe, through 
the combined use of genetic and metabolomic approaches. Among 
other compounds, we were able to identify the direct precursor of 
EGT, hercynylcysteine sulfoxide. We constructed an EGT and 
hercynylcysteine sulfoxide overproduction system utilizing three 
different overexpression strains with the inducible nmtl promoter. 
We quantified intracellular EGT content in wild type (WT) as well 
as overexpression cells by constructing a calibration curve from 
pure EGT standard injections into the LC-MS. Further, we show 
the accumulation of EGT under long-term quiescence. Finally, we 
demonstrate that the EGT pathway can also synthesize seleno- 
neine, a selenium-containing derivative of EGT, the production of 
which involves a novel intermediate compound, hercynylseleno- 
cysteine. 

Results 

EGT biosynthesis pathway 

The reported EGT biosynthetic pathways in M. smegmatis and jV. 
crassa are schematized in Figure lA. On the basis of sequence 
homology, the S. pombe locus, mugl58^ /SPEC 1604.01, was 
previously suggested to encode the main EGT biosynthetic 
enzyme [12,13]. This enzyme, with 773 amino acids, catalyzes 
triple methylation of histidine to hercynine („Ab:,A"a,Aa-trimethyl 
histidine) and subsc'qu(;nt conjugation with cysteine and oxygen to 
form hercynylcysteine sulfoxide. m!<g75iy*"/SPBCl 604.01 is a 
distant homolog of the mycobacterial EgtD and EgtB genes, 
encoding a single fusion protein. Figure IB shows the domain 
structure of this protein, according to the Conserved Domain 
Database [23], in comparison with its homologs in Schizo.scwchar- 
omyces japonicus (a relative of S. pombe), JV. crassa, and M. smegmatis. 
The exact locations and sequences of individual domains are 
shown in Table SI and the amino acid sequence alignment is 
shown in Figure SI. Due to its structure and sequence similarities 
to NcEgt-1, we presume that the S. pombe homolog also utilizes 
cysteine as a substrate, rather than using y-glutamyl-cysteine, as in 
the case of the bacterial EgtB enzyme. The bacterial pathway 
requires a subsequent removal of the glutamyl residue by another 
enzyme, EgtC [12], which has no obvious homolog in the genome 
of S. pombe. Since the mugl58^/SPRC\&M.Q\ gene encodes the 
first step of EGT synthesis, we hereafter designate it egtl'^. 

The final step of EGT biosynthesis in bacteria is represented by 
EgtE, a pyridoxal phosphate (PLP)-binding cysteine desulfurase 
that forms the end product, EGT, from hercynylcysteine sulfoxide, 
by cleaving the cysteine residue at the sulfur atom. To date, the 
EgtE homolog has not been characterized in jV. crassa. The S. pombe 
genome contains four cysteine desulfurases that could be 
homologous to EgtE: Nfsl (SPBC21D10.1 Ic), SPBC660.12c, 
SPAC11D3.10, and SPCC777.03c (Figure IC; location and 



sequences of conserved domains in Table S2; amino acid sequence 
alignment in Figure S2). Among these, SPAC11D3.10, and 
SPCC777.03c have highly similar sequences, possibly originating 
from horizontal gene transfer (Figure ID). 

To identify the EgtE homolog in S. pombe, we obtained deletion 
mutants from the Bioneer haploid deletion library [24] , cultivated 
them under nitrogen starvation (EMM2 medium lacking NH4CI, 
hereafter designated EMM2-N) for 24 hours to induce EGT 
production, and performed a metabolomic analysis. Among the 
four deletion mutants tested, ASPBC660.12c was the only one 
showing a substantial decrease in EGT and an increase in its 
precursor, hercynylcysteine sulfoxide (Figure IE; numerical results 
of all LC-MS measurements shown in Table S3). We thus 
designated the SPBC660.12c locus egt2^, as it represents the gene 
primarily responsible for the second step of EGT biosynthesis in S. 
pombe. 

Verification of egtl'^ and egt2^ by metabolomic analysis 

To verify correct assignment of the egtl''' and egt2^ genes, we 
newly constructed full deletion mutants, Aegtl and Aegt2, by 
replacing the target loci in the WT h~ 972 strain with the 
kanamycin resistance marker (kanMX). Correct integration of the 
kanMX modules into the new strains was verified by PGR. No 
difference from WT in cell size or shape was observed in these 
strains under any of the analyzed conditions. Table 1 shows the 
results of metabolome analysis of the constructed strains under 
starvation (EMM2-N or EMM2 medium with low concentration 
of glucose, 1.1 mM, hereafter designated EMM2-LG, for 
24 hours). Data regarding EGT and its precursors, starting with 
histidine, are shown. We confirmed the absence of all pathway 
intermediates in Aegtl, and the accumulation of hercynylcysteine 
sulfoxide in Aegt2. A small amount of EGT was still found in Aegt2 
under star\'ation. To determine whether any other cysteine 
desulfurase might contribute to this enzymatic reaction, we 
constructed multiple deletion mutants of SPBC660.12c {egt2^ 
with the other candidate EgtE homologs, but even in successfully 
constructed double and triple mutants, a significant amount of 
EGT stiU remained (Figure S3). This, however, is not a surprising 
result, as Seebeck [12] previously demonstrated that hercynylcys- 
teine sulfoxide could spontaneously convert into EGT in the 
presence of PLP. Furthermore, as also shown by Seebeck [12], this 
reaction could be catalyzed by an unrelated PLP-binding P-lyase 
originating from Erwinia tasmaniensis. Thus, we suspect the residual 
EGT found in Aegt2 might be a product of a hercynylcysteine 
sulfoxide reaction with PLP, possibly catalyzed by an unrelated 
PLP-binding enzyme [S. pombe genome contains at least 26 PLP- 
binding enzymes, according to PomBase [25]). 

Neither the Aegtl nor the Aegt2 strain showed any growth 
defects during cultivation in either rich (YE) or minimal (EMM2) 
culture media. Furthermore, deletion of egtT^ or egt2^ caused no 
significant perturbation to the intracellular metabolome of 
quiescent cells (Figures 2 A and 2B). Apart from the disappearance 
of EGT and its precursors, all other metabolite levels remained 
within the range of common experimental error. 

Conservation of egtl^ and egtT^ in other species 

As observed by Sc(;hc( k [12], the EGT biosynthestic pathway 
can be found in a relatively small numb(;r of eukaryotes, notably in 
the phyla Basidiomjcota (most species) and Ascomycota (most species 
in subphyla Pezizomycotina and Schizosaccfiarorrycetes). We summa- 
rized homologs of EGT biosynthesis genes in several organisms 
closely related to S. pombe (Table 2). The exact locations of the 
conserved domains in these enzymes are shown in Tables S4 and 
S5. According to sequence alignment (Figures S4 and S5), 
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Figure 1 . Characterization of EGT biosynthesis in 5. pombe. A. Previously published EGT biosynthesis pathways in M. smegmatis and N. crassa. 
The Egt-1 protein in N. crassa is a fusion enzyme that catalyzes two different reactions. B. Comparison of the conserved domain structure of S. pombe 
protein SPBC1604.01 (designated Egtl in this manuscript) with S. japonicus SJAG_00832 and EGT biosynthesis proteins in N. crassa (Egt-1) and M. 
smegmatis (EgtD and EgtB). All proteins are composed of: 1. an S-adenosyl-L-methionone (SAIVI)-dependent methyltransferase domain, including 
DUF2260, a domain of unknown function; 2. an uncharacterized DinB_2 domain, including an iron-binding motif HX3HXE; and 3. a formylglycine 
generating enzyme (FGE)-sulfatase domain. Percentage identity (% id) of amino acid sequences is indicated in comparison to the corresponding 
sequence in S. pombe. C. Comparison of conserved domain structure of M. smegmatis protein EgtE with its four putative homologs in S. pombe. All 
proteins contain a single pyridoxal phosphate (PLP)-binding cysteine desulfurase domain. The conserved catalytic residues (PLP binding sites) are 
indicated by red lines. Percentage identity (% id) of the amino acid sequences is indicated in comparison to the corresponding sequence in M. 
smegmatis. D. Phylogenetic tree visualizing the similarity of amino acid sequences of M smegmatis EgtE and its putative homologs in S. pombe. E. 
Normalized peak areas of hercynylcysteine sulfoxide and EGT obtained by metabolomic analysis of WT and deletion mutant S. pombe strains. Cells 
were nitrogen-starved prior to analysis (24 h in EIV1M2-N medium) to induce EGT synthesis. 
doi:10.1371/journal.pone.0097774.g001 



conserved domains show higher homology among species than 
inter-domain regions. Methyltransferase and FGE-sulfatase do- 
mains contain long non-homologous sequence inserts in U. majdis 
UM00197. The function of the inserts in these two domains is 
unknown, however, they are not unusual among various species 
[13]. 

Based on sequence homology to S. pomhe Egt2 (33% amino acid 
sequence identity), we propose that the NCU 11365 gene may 
encode the second enzyme for EGT biosynthesis in jV. crassa. Bello 
et al. [13] previously suggested the NCU04636 gene (19% identity 
to Egt2), which appears to correspond to nfsl'^, the mitochondrial 
cysteine desulfurase in S. pombe (63% identity). Both S. pomhe nfsl^ 
and jV. crwisa NCU04636 appear to be the closest homologs of 
mycobacterial EgtE (as also indicated in Figure ID). However, we 
did not detect any change in hercynylcysteine sulfoxide levels in 
the t^njsl deletion mutant (Figure IE). Furthermore, the njsl^ gene 
has apparent homologs in species that do not produce EGT, such 
as S. cerevisiae or humans (both genes called NFSl). 

Overexpression of egf 

To study the overproduction effect of EGT in S. pombe, we 
applied the method of Bahler et al. [26] to replace the egtl'^ native 
promoter with the nmtl promoter, which is inducible in the 
absence of thiamine [27]. We employed three versions of the nmtl 
promoter plasmid with increasing strength of expression and 
constructed three strains P81nmtl-egtl^ , P41nmtl-egtl*, and 
P3nmtl-egtl'^ , respectively. Using metabolomic analysis we con- 
firmed the accumulation of EGT and its precursors in these 



strains, and this accumulation was effectively suppressed by the 
addition of 5 |J,g/ml thiamine to the EMM2 medium (Figure 2C). 

Intracellular EGT content 

Wild type S. pombe cells contain only trace amounts of EGT 
under normal vegetative conditions (Table 1). However, EGT 
increases in quiescent cells under starvation [18,19]. We measured 
the areas of EGT peaks under vegetative, quiescent, and egtl - 
overexpressing conditions, and converted them to absolute 
concentrations using a calibration curve based upon pure EGT 
injections ranging from 1 fmol - 10 nmol (Figure S6). Intracellular 
volume was assumed to be 148.5 |a.m^ for vegetative cells [28]. For 
nitrogen- and glucose-starved cells, intracellular volumes were 
estimated as 1/3 and 2/3 of the vegetative cell volume, 
respectively. The resulting intracellular concentrations are shown 
in Table 3. 

To assess long-term variation in EGT content, we measured the 
level of intracellular EGT using metabolomic analysis during a 20- 
day time course under quiescence induced by nitrogen (EMM2-N) 
and glucose (EMM2-LG) starvation, respectively (Figure 2D). Cells 
showed time-dependent accumulation of EGT, despite being 
deprived of nutrients, suggesting that EGT might support cellular 
health under long-term quiescence. However, no loss of viability 
was observed in l^egtl mutant during 20 days of starvation 
(Figure 2E). 



Table 1. Normalized peak areas of the four compounds composing the EGT biosynthetic pathway obtained by metabolomic 
analysis of WT and newly constructed strains. 





Compound, peak m/z and 
retention time/Strain, 


Histidine 


Trimethyl-histidine 
(hercynine) 


Hercynyl-cysteine 
sulfoxide 


Ergothio-neine 


cultivation condition 


156.077 m/z @12.4 min 


198.124 m/z @10.3 min 


333.123 m/z @12.2 min 


230.096 m/z @12.6 min 




EMM2 


14.4 


1.7 


0 


0.1 


WT 972 


EMM2-N (24 h) 


2.2 


3.2 


2.2 


13.7 




EMM2-LG (24 h) 


55.4 


65.9 


3.5 


6.7 




EMM2 


10.7 


0.3 


0 


0 


Aegf) 


EMM2-N (24 h) 


2.8 


0.1 


0 


0 




EMM2-LG (24 h) 


54.1 


0 


0 


0 




EMM2 


11.4 


0.9 


3.9 


0 


Aegf2 


EMM2-N (24 h) 


1.9 


1.8 


58.2 


3.1 




EMM2-LG (24 h) 


61.1 


64.7 


441 


1.6 



Values were measured from metabolome samples of four different S. pombe strains in three different cultivation conditions, as indicated. Mass values (m/z) and LC 
retention times (min) of each peak are included for reference. 
doi:l 0.1 371/journal.pone.C097774.t001 
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Time after media change (days) 

Figure 2. Characterization of Aegt7 and Aegt2 strains. A and B. A scatter plot comparing results of metabolomic analysis of WT vs. Aegtl (A) or 
Aegt2 (B) strains under nitrogen starvation (24 h in EMM2-N medium). Each dot represents a single identified metabolite. Values on both scales 
indicate normalized peak areas of metabolite peaks in corresponding strains. Red diagonal lines indicate a 2-fold difference. C. Results of 
metabolomic analysis of WT and egtl* overexpression strains. Cells were cultivated at 26"C in the EIV1IV12 medium lacking thiamine for at least 24 h. 
Cultures indicated +Thiamine were cultivated in the presence of 5 |ig/ml thiamine for 24 h. Normalized peak areas of compounds composing the EGT 
pathway are shown. D. Time course metabolomic analysis of quiescent S. pombe cultures under nitrogen (EIV1M2-N) and glucose (EMIVI2-LG) 
starvation. Values represent means ± standard deviations of normalized peak areas of EGT in three independent cell cultures. E. Time course viability 
results of WT and Aegtl deletion mutant cultures under nitrogen (EMIV12-N) and glucose (EMIVI2-LG) starvation. 
doi:1 0.1 371 /journal.pone.0097774.g002 
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Table 2. Closest homologs of S. pombe Egtl and Egt2 proteins in selected species. 





Organism 


Closest homolog of 5. pombe Egtl 


Closest homolog of 5. pombe Egt2 


Scbizosaccharomyces japonicus 


SJAG_00832 


SJAG_03856 


Schizosaccharomyces octosporus 


SOCG_01424 


SOCG_02548 


Neurospora crassa 


NCU04343 (NcEgt-1) 


NCU11365 


Aspergillus niger 


Anl5g05880 


An02g02030 or An05g02190 


Aspergillus oryzae 


AO09001 2000265 


AO090026000291 


Ustilago maydis 


U MOO 197 


UM04128 


Candidate homologs were searched 
candidates with the best similarity 


using the on-line version of the Basic Local Alignment Search Tool (http://blast.ncbi.nlm.nih.gov) for protein sequences (blastp) and 
icores (lowest blastp E-values) were selected. 



doi:l 0.1 371 /journal.pone.0097774.t002 



Contribution of egtl^ to oxidative stress response 

As EGT is generally considered to be a physiological 
antioxidant, and the ANcEgt-1 mutant was reportedly sensitive to 
feri-butyl hydroperoxide in jV. crassa [13], we performed spot test 
experiments on EMM2 plates containing oxidants hydrogen 
peroxide and feri-butyl hydroperoxide, using the deletion and 
overexpression mutants described above (Figure S7). However, we 
did not observe any sensitivity or resistance of these strains 
compared to WT 972 strain, suggesting that egtl'^ might not be 
among the primary mechanisms that protect S. pomhe from 
exogenous peroxide. 

Biosynthesis of selenoneine 

Selenoneine (Figure 3A) is a selenium-containing derivative of 
EGT found in tuna, and possibly implicated in methylmercury 
detoxification [29,30]. As the EMM2 medium does not normally 
contain selenium, it is not surprising that no selenoneine was 
detected in previous S. pombe metabolome data sets. To test 
whether S. pombe can produce selenoneine, we first examined cell 
cultivation in a liquid EMM2 medium supplemented with various 
concentrations of Na2Se04 (Figure S8). The medium with 10 |J,M 
NajSeO,! (hereafter designated EMM2+Se) was selected for 
further experiments, as S. pombe cells exhibited quite normal 
(albeit slighdy slower) proliferation in this condition. Metabolomic 
analysis was performed on WT vegetative cells, WT nitrogen- 
starved cells (24 h in EMM2-N-I-Se), and P3nmtl-egtl'^ overexpres- 
sion mutant cells. We could clearly observe accumulation of 
selenoneine in the PSnmtl-egtl'^ strain cultivated in EMM2+Se 
medium (Figure 3B), suggesting that the overexpressed egtl^ gene 
was also responsible for selenoneine synthesis. A tiny signal of 
selenoneine (< 1 % of the EGT signal intensity) could also be 
detected in WT nitrogen-starved cells, and the signals of both 



EGT and selenoneine disappeared in the ^egtl deletion mutant 
(Figwe 3C). 

To rule out the possibility that selenoneine was produced by 
direct conversion from ergothioneine (without requiring any egtl^ 
activity), we performed two additional experiments. First, pure 
EGT was mixed with an equal concentration of Na2Se04 in vitro 
and incubated at room temperature for 24 h. No selenoneine 
signal was detected in this mixture (Figure S9). In the second 
experiment we supplemented the lS.egtl mutant with 1 mM EGT, 
which was apparently transported into the cells and produced a 
strong EGT peak. No selenoneine was detected in this case either 
(Figure 3D). Furthermore, the weak, but clearly detectable 
selenoneine signal found in WT cells did not increase as a result 
of EGT supplementation. These results suggest that egtl'^ activity is 
indispensable for selenoneine biosynthesis in S. pomhe. 

As we did not observe any signal of the presumed intermediate, 
hercynyl-selenocysteine sulfoxide, we constructed a double mutant 
of P3nmtl-egtl and Aegt2. This mutant should accumulate large 
amounts of the intermediate, assuming that the Egtl/Egt2 
pathway is used for selenoneine synthesis. Surprisingly, we found 
a strong signal of hercynylselenocysteine (not sulfoxide) in this 
mutant (Figure 3E). This signal was also found in the P3nmtl-egtl'^ 
single mutant, but further increased in the double mutant with 
lS.egt2. We thus conclude that selenoneine biosynthesis, unlike EGT 
biosynthesis, does not produce a sulfoxide as its intermediate, but 
produces hercynylselenocysteine instead. 

Since selenium naturally occurs in a very characteristic set of 
isotopes, we verified the identities of selenoneine and hercynylse- 
lenocysteine peaks by checking their isotope distribution patterns 
(Figure SIO). Finally, to check whether the EGT/selenoneine 
pathway could possibly be involved in selenium detoxification in S. 
pomhe, we performed a spot test experiment on EMM2 plates 



Table 3. Absolute intracellular EGT concentrations ( ijIVI) in S. pombe cells. 





Cell condition 


Culture medium 


Intracellular EGT ( |iM) 


WT vegetative 


EMM2 


0.3 


WT nitrogen starvation 


EMM2-N (24 h) 


157.4 


WT glucose starvation 


EMM2-LG (24 h) 


41.6 


P8lr)mt1-egtl* 


EMM2 


32.4 


P4lnrr)tl-egtl* 


EMM2 


181.2 


P3nmtl-egtl* 


EMM2 


1 606.3 



Intracellular concentrations were derived from measured normalized peak areas using a calibration curve generated by Injections of pure EGT in 10-fold dilution steps. 
The detailed calculation method is described in Figure S6. 
doi:l 0.1 371 /journal.pone.0097774.t003 
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Figure 3. Production of selenoneine in 5. pombe. A. Chemical structure of selenoneine. B. Results of nnetabolomic analysis of WT and PSnmtl- 
egtl* strains in EMM2+Se and EMiV12-N+Se nnedia. Nornnalized peak areas of EGT and selenoneine are shown. C. Extracted ion chromatogranns of EGT 
and selenoneine masses in raw LC-MS data acquired from metabolome samples of nitrogen-starved cells (24 h in EMM2-N+Se medium) of WT and 
Aegtl strains. Note that the intensity scale of the selenoneine plot is 1% relative to that of the EGT plot. D. Extracted ion chromatograms of EGT and 
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selenoneine masses in raw LC-MS data acquired from metabolome samples of WT and Aegtl strains cultivated in EMIV12+Se medium with and 
without supplementation with 1 mlVl pure EGT. Note that the intensity scale of the selenoneine plot is 0.1% relative to that of the EGT plot. E. 
Extracted ion chromatograms of six compound masses in WT, P3nmtl-egt1*, and P3nmtl-egtl* Aegt2 strains. The plots of hercynylcysteine and 
hercynylselenocysteine sulfoxide show mass values calculated from their predicted chemical formulas (C,2H2iN404S* =317.128 m/z for 
hercynylcysteine, and Ci2H2iN405Se* =381.067 m/z for hercynylselenocysteine sulfoxide, respectively). Note that the intensity scale of the plots 
in the right half of the figure is adjusted to 10% relative to the plots in the left half of the figure. 
doi:10.1371/journal.pone.0097774.g003 



containing various concentrations of Na2Se04. However, none of 
the analyzed mutants showed any growth differences compared 
with WT cells (Figure SU). 

Discussion 

In this study we applied metabolomic analysis to identify the 
egll^ and egt2^ genes composing the EGT and selenoneine 
biosynthetic pathway in S. pombe (Figure 4). These two genes have 
no homologs in S. cerevisiae, consistent with the fact that budding 
yeast do not produce EGT. In the future, a comparison between 
these two yeasts that appear similar, but are genetically rather 
distant, might provide useful clues regarding the native physiology 
of EGT or selenoneine. The presence or absence of the EGT 
pathway could be related to differences in the ecology of these two 
species. 

In an environmental stress study, transcription of the egtl'*' gene 
was up-regulated ~ 10-fold in the presence of 0.5 mM HgOj, and 
~4-fold in the presence of Cd^"*" [31]. It thus appears that egtl'*' is 
strongly stress-responsive, providing support for the hypothesized 
antioxidative role of EGT. In addition, the egtl'^ locus was 
previously named mugl58^, denoting Meiotically Up-regulated 
Gene 158, due to its up-regxilation upon entry into meiosis [32]. 
The egtf' promoter contains the TR-box motif 5'- 
TTCTTTGTTY-3', recognized by the sexual development 
transcription factor, Stell [33], and a 5'-GTAAAYA-3'-binding 
motif recognized by the forkhead transcription factor, Mei4 [34] . 
Indeed, the egtl'^ transcript was up-regulated in strains overex- 



pressing stell'^ or mei4^ [35,36]. Also, egtl^ was up-regulated in a 
strain overexpressing both basic leucine zipper (bZIP) transcrip- 
tion factors atf21^ and atf31* [36]. Furthermore, egtl^ was 
previously reported as a putative regulatory target of the 
transcription factors, Ams2, Php5, and Perl, by analysis of 
multiple genome-wide microarray datasets [37]. Ams2 is a cell 
cycle-dependent GATA factor activated during replication [38]. 
Php5 is a subunit of the CCAAT-binding complex, which 
regulates respiration, based on glucose and iron availability [39]. 
Perl is a CREB/ATF protein involved in stress response and 
sexual development [40,41]. In addition, the 5' UTR region of 
egtt contains the CuSE element sequence, 5'-DWDDHGCTGD- 
3', which is recognized by the Cufl transcription factor and 
activated by copper deficiency [42]. However, no transcriptional 
activation of egtl'^ was found under varying copper levels [43], 
suggesting this 5 ' UTR DNA fragment might not correspond to an 
actual Cufl binding site. In conclusion, it seems that egtl'^ might be 
under regulatory control of a variety of different transcription 
factors. 

In our study, the I^egtl deletion mutant showed a complete 
absence of EGT and all its precursors. In the lS.egt2 deletion 
mutant, some amount of EGT remained, consistent with results 
previously reported by Seebeck [12], that hercynylcysteine 
sulfoxide can spontaneously convert to EGT in the presence of 
PLP. Judging from the wealth of published transcriptome data, 
transcription of egt2 does not vary in response to environmental 
conditions [31], and is only mildly up-regulated (~2-fold) during 
meiosis [32], suggesting that egtl'^ represents the main regulatory 
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Figure 4. Summary of the described EGT and selenoneine biosyntKietic pathway in 5. pombe, and its transcriptional regulation. 

doi:1 0.1 371/journal.pone.0097774.g004 
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element in this patliway. Tlie Egt2 enzyme could be ubiquitously 
present in the cells, simply converting any available hercynylcys- 
teine sulfoxide into EGT. 

Overexpression strains, PSlnmtl-egtl'^, P41nmtl-egtl'^ , and 
P3nmtl-egtl^ , showed gradual accumulation of EGT in standard 
EMM2 medium without thiamine, providing an additional 
confirmation of the correct egtl* assignment. Interestingly, 
PSnmtl -egtl'*' cells contained a large amount of hercynylcysteine 
sulfoxide. A possible explanation is that the activity of Egtl in this 
strain was higher than the native acti\'ity of Egt2, resulting in 
accumulation of the metabolic intermediate between these two 
enzymes. These overexpression strains wiU be invaluable for 
verification of any proposed EGT mechanism. We also demon- 
strated that the P3nmtl-egtl'^ strain could synthesize a considerable 
amount of selenoneine when selenium was supplemented in the 
culture medium. Interestingly, selenoneine biosynthesis, unlike 
EGT biosynthesis, did not involve a formation of a sulfoxide 
intermediate, but rather involved a simple conjugate compound, 
hercynylselenocysteine. It is unknown whether egZ7''"-mediated 
synthesis of hercynylselenocysteine requires oxygen, as in the case 
of EGT synthesis. The P3nmtl-egtl strain proved to be very useful 
to elucidate the complete pathway, and we propose that this strain 
could also be employed industrially to produce EGT, as well as 
hercynylcysteine sulfoxide, hercynylselenocysteine, or selenoneine 
(the latter three compounds are not commercially available at 
present). The signal intensity of selenoneine found in WT 
nitrogen-starved cells was rather low; however, selenoneine might 
potentially have an interesting function in S. pombe. Selenoneine 
can be found in humans [44] , and it was suggested that its radical- 
scavenging activity is cxc.ii higher than that of EGT [30]. 
Importantly, we showed that selenoneine does not seem to be 
involved in detoxification of selenium in S. pombe. 

From a wider perspective, the ubiquitous presence of EGT in 
living organisms, from bacteria to humans, suggests a crucially 
important function, yet this function has not been convincingly 
demonstrated so far. We found only trace amounts (0.3 JtM) of 
EGT in vegetatively growing S. pombe cells; however, under 
starvation-induced quiescence it increased several hundred-fold. 
That, together with reported up-regulation of the egtl'^ transcript 
in the presence of II7O2 and in meiosis, implies a supportive or 
protective role of EGT for long-term hibernation of cells. 
However, the /^egtl deletion strain did not show any obvious 
defects in either sporulation or in quiescence. Furthermore, we did 
not detect any sensitivity of the l^egtl mutant to common oxidants 
H2O2 and fert-butyl hydroperoxide. As cells contain multiple 
redundant systems to deal with oxidative stress, it is possible that 
the lack of EGT could be compensated by another mechanism. 
On the other hand, the observation that the highly overexpressed 
P3nmtl-egtT^ strain, which accumulated EGT well beyond normal 
physiological levels, did not acquire any resistance to the tested 
oxidants, is intriguing. An attractive possibility exists that the true 
physiological purpos(; of EGT might lie in a yet unexplored - and 
possibly unexpected - area. We propose that genetic and 
metabolomic analyses, together with the collection of S. pombe 
strains introduced in this manuscript, may provide ideal tools to 
further investigate the in vivo role of this enigmatic compound. 

Materials and Methods 

Amino acid sequence alignment and bioinformatic 
analysis 

Amino acid sequence alignment was performed using the on- 
line version (http://www.ncbi.nlm.nih.gov/tools/cobalt/) of the 
Constraint-based Multiple Alignment Tool (COBALT) [45] with 



default parameters. A phylogenetic tree was generated from the 
on-line COBALT tool (using the Fast Minimum Evolution method 
with Max Seq Difference set to 0.9 and other parameters set to 

default values) and visualized using Analyses of Phylogenetics and 
Evolution (APE) software [46] . Percentage identity of the amino 
acid sequences was calculated using the Clustal-Omega algorithm 
[47] with default parameters. 

Strains and growth conditions 

The S. pombe strains used in this manuscript are listed in Table 4. 
The synthetic minimal medium (EMM2), rich yeast extract 
medium (YE) and sporulation-inducing medium (MEA) recipes 
were used as pubhshed previously [48] . The following variants of 
the liquid EMM2 medium were used: EMM2-N (EMM2 lacking 
NH4CI), EMM2-LG (EMM2 containing 1.1 mM - or 0.2 g/1 - 
glucose), EMM2+Se (EMM2 containing additional 10 |J,M 
Na2Se04), and EMM2-N+Se (EMM2-N containing additional 
10 |a,M Na2Se04). Cell cultures were cultivated at 26°C. 

Construction of mutants 

All DNA recombinant strains were constructed using a two-step 
PGR method. In the first step, two approximately 500-bp regions 
were amplified using genomic DNA of the WT 972 strain as a 
template, corresponding to the forward and reverse ends of the 
recombination cassette. In the second step, both modules were 
combined with the appropriate plasmid containing the kanamycin 
resistance marker (kanMXG). Transformants were selected by 
resistance to geneticin (G418) and correct integrations were 
verified by PGR. Primer sequences used for all PGR amplifications 
are included in Table S6. 

The gene disruption strains (TP1770 and TP1771) were 
constructed by replacing the target open reading frames with the 
kanamycin resistance marker. The pFA6a-kanMX6 plasmid [26] 
was used as a template for construction of replacement cassettes. 
The overexpression strains integrating the nmtl* promoter 
(TP18,57, TPhS.-^.^i, and TP1803) were constructed using the 
pFA6a-kanMX6-P81nmtl, pFA6a-kanMX6-P41nmtl, and pFA6a- 
kanMX6-P3nmtl plasmids [26] as templates. 

Hygromycin-resistant (hphMX6) and clonNAT-resistant 
(natMX6) versions of the deletion strains (TP1707, TP1732, and 
TP1813) were prepared by PGR amplification of a marker switch 
cassette and transformation of the cassette into the original strains, 
as described previously [49]. The pAG32 and pAG25 plasmids [50] 
were used as templates for the construction of the marker 
switching cassettes. Crosses were performed by sporulation on 
an MEA plate followed by tetrad dissection of the formed asci and 
selection on YE plates containing the appropriate combination of 
drugs. 

Metabolome sample preparation 

Metabolomic analysis was performed as previously described 
[17]. Briefly, cells from cultures (40 ml/ sample, 3.3x10 cells/ml 
for vegetative cells, or 10^ cells/ml for nitrogen-starved cells, 
respectively) were collected by vacuum filtration and immediately 
quenched in 25 ml of — 40°C methanol. Cells were harvested by 
centrifugation at — 20°C and constant amounts of internal 
standards (10 nmol of HEPES and PIPES) were added to each 
sample. Cells were disrupted using a Multi-Beads Shocker (Yasui 
Kikai, Osaka, Japan). Proteins were removed by filtering on an 
Amicon Ultra 10-kDa cut-off filter (MUlipore, BUlerica, USA) and 
samples were concentrated by vacuum evaporation. Finally, each 
sample was re-suspended in 40 |il of 50% acetonitrile and 1 |J,1 was 
used for LC-MS analysis. 
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Table 4. S. pombe strains used in this manuscript. 




Strsin nsinc 


G6notyp6 


Source 


972 


(WT) 


Leupold (1950) [52] 


975 


h* (WT) 




KS1366 


Astyl::ura4* ura4-D18 


Shiozaki and Russell (1995) [53] 


TP1701 


hnfsT::kanMX4 


Strains from the Bioneer haploid deletion mutant collection [24] 
were backcrossed with WT 972 to remove auxotrophic markers. 


TP1704 


ASPBC660.12c::kanMX4 




TP1705 


h" &SPACllD3.10::kanMX4 




TP 1706 


h* ASPCC777.03c::kanMX4 




TP1707 


h" ASPACUD3.10::hphMX6 


Marker switch of TPl 705 to hphMX6 


TP1732 


h" &SPBC660. 12c::natMX6 


Marker switch of TP 1704 to natMX6 


TP1733 


h* ASPBC660.12c::natMX6 


TPl 732 crossed with WT 975 


TP1736 


h" &SP8C660.12c::natMX6 Anfsl :.kanMX4 


TPl 701 crossed with TPl 733 


TP1737 


h" &SPBC660. 12c::natMX6 ASPCC777.03c::kanMX4 


TPl 706 crossed with TPl 732 


TP1739 


h* ASPBC660.12c::natMX6 ASPACl lD3.10::hphMX6 


TPl 707 crossed with TPl 733 


TP 1740 


&SPCC777.03c::kanMX4 ASPACl 1D3. 10::hphMX6 


TPl 706 crossed with TPl 707 


TPl 743 


ASPBC660.12c::natMX6 ASPACllD3.W::hphMX6 ASPCC777.03c::kanMX4 


TPl 739 crossed with TPl 740 


TP1770 


b" Aegtl::kanMX6 




TPl 771 


b" Aegt2::kanMX6 




TPl 857 


b^ egtl::P81nmtl-egtl* 


Constructed as part of this study. 


TPl 855 


b" egtl::P41nmtl-egtl* 




TPl 803 


b" egtl::P3nmtl-egtl* 




TP1813 


b" &egt2::bphMX6 


Marker switch of TPl 771 to hphMX6 


TP1814 


b* Aegt2::bphMX6 


TPl 81 3 crossed with WT 975 


TPl 879 


b" egtl::P3nmtl-egtl* Aegt2::hphMX6 


TPl 803 crossed with TPl 81 4 


doi:1 0.1 371 /journal.pone.0097774.t004 



LC-MS analysis 

LC-MS data were obtained using a Paradigm MS4 HPLC 
system (Michrom Bioresources, Auburn, USA) coupled to an LTQ 
Orbitrap mass spectrometer (Thermo Fisher Scientific, Waltham, 
USA). LC separation was performed on a ZIC-pHILIC column 
(Merck SeQuant, Umea, Sweden; 150x2.1 mm, 5 |J,m particle 
size). AcetonitrUe (A) and 10 mM ammonium carbonate bulfer, 
pH 9.3 (B) were used as mobile phase, with gradient elution from 
80% A to 20% A in 30 min at a flow rate of 100 \xl/mm. Peak 
areas of metabolites of interest were measured using the MZmine 
2.10 software [51] and normalized by the weighted contribution of 
the peak areas of the spiked internal standards (HEPES and 
PIPES). Identification of metabolites reported in this manuscript 
was based on their theoretical m/ z values and MS/MS 
fragmentation data, with the exception of EGT, the retention 
time of which was verified by analyzing a pure standard (obtained 
from Tetrahedron, Vincennes, France). The identity of selenium- 
containing compounds was verified by their isotope distribution 
patterns (Figure SIO). 

Spot test assays 

Plate media were prepared using standard EMM2 recipe [48] 
with 17 g/1 agar. Reagents were added to the media after 
autoclaving to final concentrations as indicated. Cells were 
cultivated to a concentration of 5x10^ cells/ ml and serially 
diluted in 6 steps (5-fold dilution in each step). 5 |J,1 of the diluted 
cultures was plated in each spot. Plates were incubated at 26°C for 
6 days. 



Viability measurement 

Cell viability was measured by plating approximately 300 
cells on a YE agar plate, incubating the plate at 26°C for 4-5 
days, and counting the number of colonies formed. Viability 
was calculated as the percentage of the number of formed 
colonies against the number of colonies formed at the first time 
point. 

Supporting Information 

Figure SI Amino acid sequence alignment of S. pombe 
SPBC1604.01, S.japonicus SJAG_00832, N. crassa Egt-1, 
M. smegmatis EgtD, and M. smegmatis EgtB proteins. 

Alignment was generated using the COBALT algorithm. Con- 
served domains are indicated according to their location in S. 
pombe SPBC1604.01. 
(TIF) 

Figure S2 Amino acid sequence alignment of M. 
smegmatis EgtE protein and its four putative homologs 
in S. pombe. Alignment was generated using the COBALT 
algorithm. The conserved catalytic residue (PEP binding site) is 
indicated by a red arrow. 
(TIF) 

Figure S3 Normalized peak areas of EGT and its 
precursors obtained by metabolomic analysis of WT, 
the \SPBC660.12c single deletion mutant, and multiple 
deletion mutants with other putative homologs of 
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mycobacterial EgtE. Cells were nitrogen-starved prior to 
analysis (24 h in EMM2-N medium) to induce EGT synthesis. 
(TIF) 

Figure S4 Amino acid sequence alignment of S. pombe 
SPBG1604.01 (Egtl) protein and its closest homologs in 
selected species. Alignment was generated using the COBALT 
algorithm. Conserved domains are indicated according to their 
location in S. pombe Egtl. 
(TIF) 

Figure S5 Amino acid sequence alignment of S. pombe 
SPBC660.12c (Egt2) protein and its closest homologs in 
other species. Alignment was generated using the COBALT 
algorithm. Conserved domains are indicated according to their 
location in S. pombe Egt2. The conserved catalytic residue (PLP 
binding site) is indicated by a red arrow. 
(TIF) 

Figure S6 Absolute quantification of EGT content in 
cells. A calibration curve was constructed by performing LC-MS 
injections of pure ergothioneine in 10-fold dilution steps, 
containing a constant amount of HEPES and PIPES standards 
(250 pmol each) for normalization (upper panel). Normalized peak 
areas were plotted against injected amounts and a regression curve 
was generated using Microsoft Excel. The formula to calculate 
absolute amount (fmol) from normalized peak area was derived 
from the regression curve formula (middle panel). Normalized 
peak areas of EGT were converted into absolute intracellular 
concentrations using estimated average cellular volumes (bottom 
panel). 
(TIF) 

Figure S7 Spot test results on hydrogen peroxide and 
tert-butyl hydroperoxide agar plates. WT, deletion, and 
overexpression strains described in this manuscript were serially 
diluted and grown on EMM2 plates supplemented with increasing 
concentrations of oxidants hydrogen peroxide (H2O2) and krt- 
butyl hydroperoxide (/-BOOH). The stress-sensitive Aj-^i strain 
was used as a positive control. 
(TIF) 

Figure S8 Cell number increase in liquid EMM2 
medium supplemented with Na2Se04. Relative cell number 
increase in 24 h was measured in liquid EMM2 medium 
supplemented with increasing concentrations of Na2Se04. Cell 
cultures were incubated at 26°C. 
(TIF) 

Figure S9 Analysis of a mixture of EGT and selenium in 
vitro. Extracted ion chromatograms of EGT and selenoneine 

masses are shown for 1 mM EGT, 1 mM Na2Se04, and mixture 
of both, incubated at room temperature for 24 h. Note that the 
intensity scale of the selenoneine plot is 0.1% relative to that of the 
EGT plot. 
(TIF) 

Figure SIO Verification of the identity of selenoneine 
and hercynylselenocysteine by their isotopic patterns. 
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